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1 Introduction 



A novel approach to the massive Yang-Mihs gauge theory has been proposed 
[U, where the divergences are consistently removed in the loop expansion. 
The removal strategy follows close the method used recently for the nonlinear 
sigma model [2]. It consists in the subtraction of the pure pole parts of 
the properly normalized one-particle-irreducible amplitudes regularized in 
D dimensions. The mathematical tools are standard, but the proof that the 
technique is consistent is rather involved [2], [3], [Ij. 

Although the subtraction method is consistent with the Slavnov Taylor 
Identities, locality and a new, ad hoc derived. Local Functional Equation, 
the perturbative series seems to be inadequate for high energy processes, 
thus casting some doubts on the validity of unitarity (although SS^ = 1 or- 
der by order in the loop expansion). It has been recently suggested [5] that 
the cause of this is due to some singularities (phase transitions) present in 
the parameters space (/3 := -^,m?). According to this scenario one can ap- 
proach the theory with the usual perturbative loop expansion for low-energy 
processes, while the high energy processes are described by the massless 
Yang-Mills theory with no remnants of the longitudinal polarizations. The 
transition between the two regimes may be studied by the lattice simulation. 
This is attempted in the present paper, where we try to study the model in 
the non perturbative regime. 

The present work was motivated by the arguments just outlined, but 
its contents and results are valid independently from them. In fact it opens 
new perspectives on the lattice gauge theories. 

An intensively studied lattice gauge model [6]- [13] turns out to be the 
perfect tool for the simulation of the massive Yang-Mills (i.e. mass a la 
Stiickelberg) . We confirm the existence of a Transition Line (TL) which 
separates a confined phase from one with physical vector boson states. The 
phase LT has an end-point around ( (3 ~ 2.2, ~ 0.381): for smaller 
j3 there is a smooth transition (crossover) from one phase to the other, 
while for larger (3 there are numerical indications of singularities in the 
derivatives with respect of m? and of /3 of the energy and of the order 
parameter (the derivative of the free energy). The deconfined phase is 
studied by using the correlators of gauge invariant fields. This allows a full 
gauge invariant approach to the model. We give numerical evidence of the 
existence of iso- vector modes for the spin one (no spin zero is present). For 
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the isoscalar fields there is a faint, but persistent, signal of an energy gap 
both for spin one and zero. Far from the LT these excitations in the iso- 
scalar channels are compatible with the threshold of two iso- vector spin one 
modes. However near the TL the energy gap in the isoscalar channels is 
lower than the threshold, thus suggesting the existence of both spin one and 
spin zero bound states. This effect happens in a band attached to the TL: 
for large /3 (i.e. higher than the end point value) the band is very narrow, 
while in the crossover region (low /3) the onset of bound states is smooth and 
on a wider region. The tantalizing question is whether this interesting region 
of the phase space will ever be reached by experiments and the presence of 
bound states confirmed. 

2 The Lattice Model 

The field theory (for the SU{2) group) in the continuum is [1] 

SyM = ^ ( - \Ga,AA]G^/[A] + ^{Aa, - Fa,f) , (1) 

where in terms of the Pauli matrices 

= y^a^, = y^aM ■= i^d^^l (2) 

0(x) is an element of the SU{2) group, parameterized by four real fields 

Q = (j)Q+ iTa(t)a-, =^ (t)l + (j? = 1. (3) 



We have 



Faf, = 2{(f)od^(pa - df,(j)Q(l)a + eabcdfiMc) ■ (4) 



The action in eq.([T]) is invariant under ql^x) G SU{2)l local-left and qr E 
SU{2)ji global-right transformations 



n'{x) = gL{x)n{x) ( c)'( \ - Of \ ^ 



The theory is not renormalizable due to the non-polynomial dependence on 
$ explicit in eq. (H]). We refer to the Refs. [1] and [2] on the new strategy 
suggested for the consistent subtraction of all the ultraviolet divergences. 
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The lattice model is constructed by assuming a nearest neighbor interac- 
tion and by requiring a naive mapping into the action ([T]) in the limit of zero 
lattice spacing. It is very important to construct a model invariant under 
the discretized versions of eqs. ([5]). The link variable is taken to be 

U{x,fj.) ~ exp(— mA^(x)). (6) 

Both U{x,fi) and the site variable i}{x) are elements of SU{2). Thus the 
action is (/? = p-) 

Se =^ D\e^Tr{l - Uo)+^M^a^y{e^Tr^l - n{x)^U{x, f^t)n{x + ^)}, (7) 
n xfi 

where the sum over the plaquette is the Wilson action [11] and the mass 
term has the (Euclidean) continuum limit 

/3 



Xfl 

j d'^xTr |(^^ - ind^n^^y (8) 



Xfl 



Se yte^Tr{Un) - ^m^^Re^Tr^n{x)^U{x, n)n{x + fi)^ (9) 



In the simulation the Tr{l} is omitted. Thus the action becomes 

From now on the dimensionless parameters are /3 and tv?. We work in 
Z? = 4, however the symbol D is kept in some equations. In the paper we 
will consider also the model with — )■ — m^. 



3 Simulation 

The partition function is obtained by summing over all configurations given 
by the link variables and the gauge field 

Z[/3,m2,iV] = ^ e-^^, (10) 
{u,n} 

where is the number of sites. 

In principle the integration over il(x) is redundant, since by a change of 
variables {Un{x,fj,) := ^l{x)^U{x, fM)i}{x + /.t)) we can factor out the volume 
of the group. Z[/3,m^,iV] becomes 

[ E ] E E + ^"^'^^ E Tr{U{x, ^)}) . (11) 

{n} {U} □ xij. 
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In eq. ([TT]) the integration over Q has disappeared; consequently Q in eq. 
(jlUp does not describe any degree of freedom. In that respect we are at 
variance with other approaches to the same action ([7|) as in [6]- |13] . where 
the field is thought of as a Higgs field with frozen length. In eq. (jlOp we 
force the integration over the gauge orbit Uq by means of the explicit sum 
over In doing this we gain an interesting theoretical setup of the model; 
in practice, our formalism is fully gauge invariant (Section [5|). Moreover 
by forcing the integration over the gauge orbit Un we get results which are 
less noisy than those obtained by using only the integration over the link 
variables in (fTTj) . 

4 Functionals and Order Parameter 

In this model we can study the energy-per-site functional 

= ^{'neJ2Tr{Un} + m^J2^r{n^{x)U{x,fi)n{x + fi)}y (12) 

□ XfJ. 

where ( ) denotes the mean value by the Boltzmann weight of eq. (jlOp . 
Moreover we introduce the magnetization, i.e. the response to the applied 



^ = 9]^(^" E Tr{n\x)Uix, ti)n{x + /x)}) . (13) 



DNB dm? 2ND 



Then we have the plaquette energy 

\ 2 r „ 1 

(14) 

There are some simple properties that will be of some help in the sequel. 
Under the mapping 



U{x,ti) ^ -Uix.fi) (15) 

the Wilson action is invariant while the mass part changes sign. The measure 
of the group integration is invariant, then we have from eqs. (jlOp . (jl2p and 

m 

Z[/3,-m^,N] = Z[l3,n?,N] 

£[/3,-m2,iV] = - e:[/3,m2,iV]. (16) 
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We argue that C is a valid order parameter. This will be shown in the next 
Sections. We briefly recollect the main points. Its suscpetibility (d^/dm?) 
has a cusp-like behavior in on the TL for /3 greater than the end point 
value. It is odd under the change of sign of m^. Moreover the implementa- 
tion of the global SU(2)ji symmetry is drastically different in the far away 
regions where £ ~ and Cit 1. In the first region SU{2)ji charged fields are 
screened or confined, while in the second deconfined modes are present. 



5 The Vector Meson Fields 

Our approach allows the presence of SU{2)l gauge invariant fields. Let us 
consider 

C{x, fi) := il.\x)U{x, n)il.{x + fJ.) = Cq{x, /i) + ITaCaix, /i), (17) 

which, according to eqs.([5]), are invariant under local SU{2)l transforma- 
tions. Cq{x, ji) is the mass term density in the action ([9]) and it is a SU (2)/j - 
scalar (isoscalar), while Ca(x, ^) are vectors under the same group of trans- 
formations (isovectors). Since C(x,/i) G SU{2), we get that all fields are 
real and constrained by 

Coix,^l)^ + Y,Caix,^J.)^ = l. (18) 

a 

Moreover we expect the vacuum to be invariant under SU{2)ji global trans- 
formations ([5]) and therefore 

(Ca(x,^))=0, a = 1,2,3, V(x,//). (19) 

The order parameter ()13p 

^=^E(^o(^'^)) (20) 

is the conjugate of the mass parameter m^. According to eq. ()16p we 
expect that 

lim C = 0, lim C = 1. (21) 

Beside the order parameter, it is important to study the following connected 
correlators. They will provide the essential characterization of the decon- 
fined phase of the system. 



c 
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Cob,^luix,y) ■■= {Co{x,fi)Cb{y,i'] 



c 



Coo 

y):=(Co(x,/i)Co(y,z.))^. (22) 

In order to investigate the transition between phases we consider also the 
susceptibihty 

^^=7^EE(^o(x,M)Co(y, 



= :^EE ((^o(^"")^o(2/,^)) - {Co{x,f,))(Co{y,u))). (23) 

It should be noticed that the mean square error of C is related to its deriva- 
tive 



This relation is very important for numerical simulations. If the derivative 
of £ had a finite limit for N ^ 00, the standard deviation would have a 
l/\/iV behavior. If instead the derivative diverges then the standard error 
might not have a decreasing behavior by increasing the lattice size A^. If 
this is the case, then the calculation of the derivative by using the heat bath 
yields a noisy signal. The noise might not decrease by increasing the lattice 
size, as expected in the normal case. 

5.1 The SU{2) Right Symmetry 

If the SU{2)fj symmetry is unitarily implemented then we expect 

Cab,tiu{x,y) = 0, if a/ 6, a,6 = 1,2,3 

CobM^^y) = 0- (25) 

The energy gap in the correlator in 6*00,^1/(2;, 2/) might set on above the 
two-particle threshold. However there is an interesting possibility that the 
gap (both for spin one and spin zero) shows up below this threshold, thus 
suggesting the existence of bound states. 

5.2 The Continuum Limit of C 

By a similar argument used in Section [2] we study the continuum limit of 
C{x, fi). We have 

C{x,n) = n\x){l - iaA^{x)){n{x) + adf,n) + 0{a^). (26) 
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Thus for Ci , C2 , C3 one gets 

iTaCa{x,^i) = -iaQ^ (^Af,{x) - md^n^^n + Oia"^). (27) 

While for Cq we can use the result of Section [2l eqs. d?]) and ([8]) 

2 

Co{x,fi) = 1 - ^Tr [{A^ - md^n^f] + 0{a'). (28) 

Notice that the dominant terms in eqs. (p7|l and (f28|) are 5C/ (2) local-left- 
invariant. 

The continuum limit of Co(x,//) and the expression of the order param- 
eter in eq. (|13p suggests that in the region C ~ where confinement (or 
screening) occurs there is condensation of vector meson pairs, while in the 
deconfined region, where C ~ 1, vector mesons do not condense and are 
realized as particles. 

6 Note on Symmetry Breaking 

The symmetry of the model and of the partition function is rather interest- 
ing. The SU (2) L left transformations ^ correspond to the local symmetries 
of the action, while the SU{2)r transformations ([5]) can only be global sym- 
metries, due to the fact that the mass term in Se (eq. ([7])) breaks the local 
SU(2)ji symmetry. For decreasing mass parameter m? we expect the onset 
of a local SU{2)ji symmetry. Then the fields C{x,fi), by construction p!7|) . 
transform according to a SU{2) (X" SU{2) ~ 0(4) group of transformations 

Cix, fiY = gR{x)C{x, ij)gR{x + (29) 

This fact has far reaching consequences. That is, in the limit of zero mass, 
only closed loops of C{x,fi) fields have non zero expectation value [15j . In 
particular all the correlators in eqs. (|22p are zero unless y = x + ^ and 
y + V = X, i.e. v = —fj. and C{y,v) = C(x,/u)^. But then the 0(4) on the 
set {Co{x),Ca{x)} imposes 

CoQ,f,v{x,y) ~ Cii^^y{x,y) = C22,^iv{x,y) = C^:i^^y{x,y). (30) 

The numerical simulations show that the onset of a local SU{2)ii is very 
rapid when one crosses the line of PT. It becomes smooth for small /3, after 
the end point. 

When the mass parameter becomes large the 0(4) symmetry will be lost 
and only the global SU (2)^^ invariance will survive and therefore 6*00,^1/(2;, y) 
will be substantially different from the SU{2)ji - vector components. 
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7 Survey 



We have performed standard Monte Carlo simulations for the model based 
on the action ([9]). Heat bath has been used for the updating. Normally 
we have saved a configuration every fifteen updatings for a total of 10,000 
measures. 

We considered cubic 4-dimensional lattices with periodic boundary con- 
ditions of different sizes: 4^, 6^, 8^, 12^, 16^, 24^. We have chosen the size on 
the basis of the precision needed. Typically a 6^ lattice size provides sensible 
results if {13, m?) is far away from the TL. 

We have built a bird's-eye view of the region /3 G [1,4], m? € [0, 8] of some 
global quantities of the system. The quantities studied are those described 
in Section m i.e. the energy per site E (fT2]) . the order parameter C (fT3]) and 
their derivatives with respect to m? and p as in eq. ()23p 

d ^ 



dm? 

d 
dm? 
= D (t 



E 




Co(x, ^) - {-^ ^ Co(a;, ^) 




-/5 



d 



1 



Se — { 



Se) 



Co{x, /i) - {-j= Co(x, fi)) 



XfJ, 



Notice that 
d ^ 



1 d 



Se — { 



Se) 



(31) 



9/3 ND dp 



1 

D 



iy^^^o{x,fj^)) 

^ xn 




PVN 



Se) 



(32) 



A survey on the parameters space shows a clear phase change across the line 
represented in Fig. [TJ In particular the line represents the loci where the 
dependence of the order parameter € from has a marked inflection. The 
line is stable under the change of the size, for instance from 6^ through 24"^. 
A throughout study has shown that both energy and order parameter are 
continuous everywhere including on the TL. Fig. [2]describes the dependence 
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Points of inflection of C as a function of (transition line). Ensemble 10' 



E 




1 1.5 2 2.5 3 

(3 (arrow at |3=2.2) 

Figure 1: Phase diagram 

on m? of C and its derivative. All the first derivatives have a cusp behavior 
for /? > 2.2 while they are smooth for (3 < 2.2. Fig. [3] exemplifies the 
situation for /3 = 1.5 and (3 = 3. There is some evidence for an end point at 
/3 ~ 2.2, ~ 0.381, linked to the crossover point evidenced by early works 
on SU(2)-QCD simulations [IS]. 



8 Numerical Results 



The numerical analysis of Sect. [7] confirms the results obtained in previous 
works about the TL, with some minor discrepancies, as the position of the 
end point. On the vexata quaestio, concerning the order of the PT across 
the line and for /3 > 2.2, our numerical evidence is not very conclusive, 
although we would be more in favor of a second order type. The present 
Section is devoted to the new and surprising results. They show that the 
model is indeed a faithful simulation of the massive Yang-Mills gauge theory 
and moreover that unexpected and non trivial features can be obtained in 
a region of the phase diagram unaccessible by perturbation theory. 
We study the operators (zero three-momentum) 



Ca,f,{t) := ^L= V'Ca(f,X4,/i) 



0,1,2,3, fi,u = 1,2,3,4. (33) 



X4,=t 
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C. p=2.2. Ensemble 10^ 



dC/dm^. (3=2-2. Ensemble lo'' 




"T 1 1 1 r 

+ + 



+ 



Size 12^ + 
_ Size_ 61._._x_ 



J I I L_ 



(arrow at 0.381) 



1.37 0.375 0.38 0.385 0.39 0.395 0.4 
™2 



Figure 2: C and ^ at /3 = 2.2 for size 6^ and 12^. 



In particular we consider the two-point correlators 

Cab,iiu{t) ■= (^Ca,ij.{t + to)Cb,u{to) 

Numerical simulations support the selection rules 



C 



Cab,fii/{t) 



0, a, 6 = 1,2, 3 



(34) 



(35) 



imposed by the global SU{2)fi invariance. Moreover, according to eqs. 
dm), (IH]) and ([301) we get the limit values 



lim Caa,^.(0) = 0.23V > 

m2=0 



a = 0,1,2,3. 



(36) 



while for t > in the same limit the two-point functions are expected to 
vanish according to the discussion of Section [6l 
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dC/dm^ P=1.5 Ensemble 10"^ 



dC/dm^ P=3 Ensemble 10^ and 10^^ 




8 
7 
6 
5 
4 
3 
2 
1 


0.6 0.7 0.8 0.9 1 1.1 1.2 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



' X 


1 1 1 1 . 




' : riof ^ ! " 
xJ 24'..l0^....x....: 




-r 

x + 


$ 

1 


X 

+ + + 

1 1 1 1 



m (arrow 0.86) 



m (arrow at 0.235) 



Figure 3: at /3 = 1.5 and /3 = 3 for different lattice sizes. 

The spin analysis is done by decomposing the correlators into a spin one 
and spin zero parts (dots stand for pairs of iso-indexes 00 or 11,22,33 ) 

C...,f,^{t) = - ^4m<54^) + S„Xt)6^f,6i^. (37) 

We fit the amplitudes by a single exponential form 

F{t) = a + be-^^. (38) 

A more complex analysis is not at reach with the data at hand. However 
the form turns out to be sufficient for most of the cases that have been 
considered. The energy gap A is obtained from a fit on the correlator (|34p 
evaluated on 10^ configurations. Typical results are shown in Figs. |l]and[5j 

The energy gaps have been evaluated for several values of (/3, m?). Figs. 
[6] and [7] represent the energy gaps as function of . Several features are 
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Cqo- 

0.35 



=3 m^=0.25 Size 1 2'^Ensemble lo"^ 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




i 



Xv*°'f 1 .94 +/- 0.08 
As<°'=1.7 +/- 0.06 
Vector fit 
Scalar fit 



0.5 1 1.5 2 2.5 
Time 



3 



3 3.5 4 



1. p=3 m^=0.25 Size 12'^Ensemble 10'* 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




~i 1 1 1 1 1 1 r 



Av*^ '=2.338 +/- 0.14 I ^ 

Vector fit 

Scalar: no fit : x 



0.5 1 1.5 2 2.5 3 3.5 4 
Time 



Figure 4: Correlators at /3 = 3 and = 0.25: Cqo for the isoscalar part 
and Cii for the isovector according to eq. ([57]) . TL is at m? = 0.235. 
The numerical values are the energy gaps for spin 1 (Ay\ ^y^) ^-^d spin 

present in all the cases we have considered, i) In the deconfined region and 
far from the TL the isovector correlator is due to a spin one mode with 
energy gap A ~ \m\, while the isoscalar correlator has both spin one and 
spin zero energy gaps, consistent with a two-vector-meson threshold, ii) 
Near the TL the isoscalar gaps become smaller than the threshold, thus 
suggesting the presence of bound states, iii) Across the TL there is a rapid 
increase of the gaps: within the errors all correlators vanish for t > and 
a 0(4) symmetry is restored. The change of phase is much more rapid for 
/? = 3 than for /? = 1.5 (the change of scale of in Figs. [6] and [7] should 
be properly accounted for). 

9 Conclusions 

We have investigated the deconfined phase of a massive Yang-Mills model 
by using a set of gauge invariant fields. We give evidence of a TL in the pa- 
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Cqo- P=1 .5 m =1 .0 Size 1 'T 



(0)_ 



: 2.486 +/- 0.05 
=2.609 +/- 0.09 
Vector fit 
Scalar fit 



0.5 1 1.5 2 2.5 3 3.5 
Time 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




Ci 1 . p=1 .5 m =1 .0 Size 1 1 



Ay*'' '=1.83 +/- 0.03 I ^ 

\t Scalar: no fit ^— x- 

Vector fit 



^ 



0.5 1 1.5 2 2.5 3 3.5 
Time 



Figure 5: Correlators at /3 = 1.5 and rr? = 1.0 (TL is at = 0.85). See 
caption of FigUl 

rameters space [P^m^]. An order parameter C is introduced (the response 
to the parameter) which is ~ 1 in the deconfined region (large m^), 
while ~ below the TL (small tv?). The vanishing of the order parameter 
corresponds to the condensation of pairs of vector mesons. Far from this 
line the spectrum consists of an isovector spin one meson and of two-particle 
states in the isoscalar spin one and spin zero channels. Moreover there is 
some evidence of bound states near the TL in the isoscalar channels for both 
spin states. The presence of a discontinuity line confirms the conjecture on 
the existence of two regimes: a low energy where the loop expansion is valid 
and an extreme region where massless Yang- Mills works. The last point is a 
very important step forward for the understanding of high energy processes. 
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Isoscalar, p=1 .5, Size = 12^^, Ens. 10"^ 




(arrow at m^=0.85) 



Isovector, p=1.5, Size = 12^, Ens. 10"* 




Figure 6: On the y-axis the energy gaps of isoscalar spin one and spin zero 
(left) and isovector spin one (right) as function of for fixed /3 = 1.5 (TL 
at rrP = 0.85). The energy gaps have been extracted from fits as in Figs. S] 
and [5] by using the expression in eq. ()38p . The line is Y^jm^j and it is very 
close to the spin 1 isovector energy gap in the region away from the phase 
transition. 
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